DISCRETE HAMILTON PONTRYAGIN MECHANICS AND 
GENERATING FUNCTIONS ON LIE GROUPOIDS 

Ari Stern 

We present a discrete analog of the recently introduced Hamilton- 
Pontryagin variational principle in Lagrangian mechanics. This unifies 
two, previously disparate approaches to discrete Lagrangian mechan- 
ics: either using the discrete Lagrangian to define a finite version of 
Hamilton's action principle, or treating it as a symplectic generating 
function. This is demonstrated for a discrete Lagrangian defined on an 
arbitrary Lie groupoid; the often encountered special case of the pair 
groupoid (or Cartesian square) is also given as a worked example. 

1. Introduction 

In a recent paper, Yoshimura and Marsden (2006) introduced the Hamilton- 
Pontryagin variational principle for Lagrangian mechanics. Given a smooth 
configuration manifold Q and a Lagrangian L : TQ M, this principle 
defines an action for paths in the so-called Pontryagin bundle TQ ®T*Q, 
whose elements are written {q,v,p). Critical paths for this action, in addi- 
tion to satisfying the usual Euler-Lagrange equations p = dL{q,v)/dq, also 
satisfy the the Legendre transform p = dL(q,v)/dv and the second-order 
curve condition q = v. This can be seen as a unification of two equivalent, 
but previously disparate, approaches to Lagrangian mechanics, where the 
Lagrangian is used either (a) to define Hamilton's action functional, study- 
ing its critical paths; or (b) to define the Legendre transform, using it to 
pull back the canonical symplectic structure from the cotangent bundle to 
the tangent bundle. 

There is a similar dilemma in the approach to discrete Lagrangian 
mechanics. Given a discrete Lagrangian L/j : Q x Q M (or more generally, 
L/^: G ^ M, where G ^ Q is a Lie groupoid), one can either (a) use to 
define a discrete version of Hamilton's action principle, or (b) treat as a 
symplectic generating function. While some progress has been made towards 
combining these approaches into a discrete Hamilton-Pontryagin principle 
for certain important special cases — namely, the pair groupoid Q x Q when 
Q is either a vector space (Kharevych et al., 2006; Leok and Ohsawa, 2008) 
or a Lie group (Bou-Rabee and Marsden, 2009) — this problem has not yet 
been resolved for general configuration manifolds Q, nor for the even more 
general case of a Lie groupoid G ^ Q. 
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In this paper, we present a discrete Hamilton-Pontryagin principle, which 
is shown to unify these two disparate approaches (variational principles vs. 
generating functions) for discrete Lagrangian mechanics on arbitrary Lie 
groupoids. We begin, in section 2, by giving a brief review of Lagrangian 
mechanics, including the continuous Hamilton-Pontryagin principle, as well 
as summarizing the existing frameworks for discrete Lagrangian mechanics 
on Q X Q and on Lie groupoids G ^ Q. Next, in section 3, we introduce 
the discrete Hamilton-Pontryagin principle, which is defined with respect 
to paths in the cotangent groupoid T*G ^ A*G beginning at the zero 
section; the approach is related to that used by Milinkovic (1999) in study- 
ing the Morse homology of generating functions. This variational principle 
and its solutions, which imply those of the previous approaches to discrete 
Lagrangian mechanics, are derived first in cotangent bundle coordinates and 
then given intrinsically. Finally, in section 4, wc work out the special case 
of the pair groupoid Q x Q, in particular showing how this corresponds to 
the formulation of Leok and Ohsawa (2008) when Q is a vector space. 

2. Lagrangian Mechanics, Continuous and Discrete 

2.1. Lagrangian and Hamilton-Pontryagin mechanics. Let Q be a 

smooth configuration manifold, and L: TQ — > R be a Lagrangian on its 
tangent bundle. There are two theoretically equivalent, but conceptually 
distinct, approaches to the Lagrangian mechanics of this system. The 
first, which we will call the variational approach, is to study critical paths 
q: [a,b] ^ Q for the action functional 



Such a path is critical if and only if it satisfies the Euler-Lagrange equations 

d_dL _ dL 
dt dq dq ' 

The second approach, which we will call the symplectic approach, is to use 

the Legendre transform ¥L: TQ — > T*Q, which takes {q,t)) f-^ {q,dL/dv), 
to pull back the canonical symplectic structure from T*Q to TQ. If oj 
is the canonical symplectic 2- form on T*Q, then one can define a 2-forni 
cjl = (FL)*a; on TQ. (The form is symplectic when L is a hyperregular 
Lagrangian, and presymplectic more generally.) Next, one defines the energy 
function E : TQ — ^ M given by 

E (q^ v) = ¥L {q,v) ■ V — L {q, v) . 

Finally, then, one looks for vector fields X e X {TQ) that satisfy 
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which is essentiaUy the tangent bundle version of Hamilton's equations on 
T*Q. (For further background, see Marsden and Ratiu, 1999, Chapter 7.) 

Yoshimura and Marsden (2006) showed that these two approaches can 
be unified through an expanded variational principle, which they call the 
Hamilton-Pontryagin principle. Given a path {q, v,p) : [a, b] — > TQ © T*Q, 
the Hamilton-Pontryagin action is given by 

S {q, v,p) = f [L {q{t), v{t)) + p{t) ■ {q{t) - vit))] dt. 
J a 

This is essentially the usual action functional — except, rather than sim- 
ply prescribing the second-order curve constraint q = v, one treats q and 
V as independent variables and then uses p as a Lagrange multiplier to 
enforce this constraint. Varying over paths with prescribed endpoints, so 
that 6q{a) = and 6q{b) = 0, the variation of the action is 



dS {q,v,p) ■ {Sq,Sv,Sp) 



J a 

J a 



— ■ 5q + — ■ 5v + p ■ (6q - 6v) + 6p ■ (q - v) 
oq ov 

dL \ ^ (dL \ , ; 

~Q^~Pj ( ^Pj ■ov + dp-{q-v) 



dt 



dt. 



Therefore, {q,v,p) is a critical path if and only if it solves the so-called 
implicit Euler-Lagrange equations 

dL dL 
oq ov 

which combine the Euler-Lagrange equations, the Legendre transform, and 
the second-order curve condition into a single set of equations. It should 
be noted that this principle is especially useful for studying constrained 
and other degenerate systems, and is closely connected with the generalized 
Legendre transform of Tulczyjew (1974, 1977), the generalized Hamiltonian 
dynamics formalism of Skinner and Rusk (1983), and the Dirac structures 
of Courant (1990). 

2.2. Discrete Lagrangian mechanics. The idea of discrete Lagrangian 
mechanics was put forward in seminal papers by Suris (1990) and Moser and 

Veselov (1991), among others, and a general theory was developed over the 
subsequent decade. (See Marsden and West, 2001, for a comprehensive sur- 
vey.) This work was motivated by the need to develop structure-preserving 
(e.g., symplectic) numerical integrators for Lagrangian mechanical systems 
on general configuration manifolds; the methods developed using this dis- 
crete Lagrangian framework are called variational integrators. 

For most of the work in this field, the starting point is to replace the 
Lagrangian L : TQ — >■ M by a discrete Lagrangian Lh: Q x Q E., which 
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approximates the contribution to the action integral, 

Lh{qo,qi)^[ L{q{t),q{t))dt, 
Jto 

for a time step of size h = ti — to. As with continuous Lagrangian mechanics, 
there are two typical ways to proceed, following either the variational or the 
symplectic point of view. 

The variational approach to discrete Lagrangian mechanics is as follows. 
Suppose that we specify a sequence of time steps a = to < h < ■ ■ ■ < = b, 
with equal step size h = tn+i — tn for n = 0,...,A'^ — 1. We then define a 
discrete path to be a sequence of configuration points qo, . . . ,qN G Q] this 
can be thought of as approximating a continuous path q: [a,b] — )■ Q, with 
In ~ Q{tn)- Given the discrete Lagrangian Lh'. Q x Q ^ the discrete 
action sum is defined to be 

JV-l -6 
Sh {qo, ...,qN)= {qn, qn+i) ^ j L {q{t), q{t)) dt. 

Next, taking fixcd-cndpoint variations of the discrete path, so that 5qo = 
and 6qM = 0, it follows that 

dSh (go, ■■■,qN) ■ {Sqo, ^qn) 

N-l 

= XI [^0'^'^ + ^1-^'* {qn-l,qn)] ■ Sqn- 

n=l 

Therefore, a discrete path is critical if and only if it satisfies the discrete 
Euler-Lagrange equations 

doLh {qn,qn+i) + diLh {qn-i,qn) = 0, n = l,...,N -1. 

This implicitly defines a two-step numerical integrator on Q x Q, which 
(given suitable assumptions of nondegeneracy) maps {qn-i,qn) '-^ {qn, qn+i)- 
On the other hand, the symplectic approach to discrete Lagrangian 
mechanics is to view : Q x Q — ^ M as the generating function for a sym- 
plectic map on T*Q. To do this, one defines the discrete Legendre transforms 
¥^Lh -.QxQ^ T*Q by 

^~Lh (go, gi) = -doLh (go, gi) , F+L^ (go, gi) = diLh (go, gi) • 
Therefore, we can implicitly define a map on T*Q by 

po = ¥~Lh{qo,qi) , pi = F+L^ (go, gi) , 
and if F~L/j is invertible, this defines one step of the symplectic integrator 

¥+Lh o (F-Lft)"' : T*Q ^ T*Q, (go,po) ^ (gi,Pi) • 

More precisely, the discrete Legendre transforms define the Lagrangian sub- 
manifold of {T*Q, —oj) X {T*Q,uj) generated by L^. Prom this perspective. 
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if the invertibility condition holds, then this submanifold is the graph of a 
symplectic map on T*Q (see Weinstein, 1971, 1979). 

Note that, if we perform the composition in the opposite order, we get 
the previously derived two-step method, 

{¥-Lh) o F+Lft : Q xQ xQ, g„) ^ g„+i) , 

where the discrete Euler-Lagrange equations follow automatically from the 
fact that ¥~Lh {qn, Qn+i) = F+L/^ {qn-i, (In)- 

2.3. Discrete Lagrangian mechanics and Lie groupoids. Weinstein 

(1996) observed that both approaches in the previous section can be gener- 
alized using Lie groupoids. Let G ^ Q be a given Lie groupoid, and define 
a discrete Lagrangian L/ji G ^ M. The earlier formulations then coincide 
with the special case G = Q xQ, which is called the pair groupoid. This per- 
spective has continued to bear fruit in recent years, being further developed 
by Marrero et al. (2006) and extended to discrete nonholonomic Lagrangian 
mechanics by Iglesias et al. (2008). 

In the variational approach, one begins by taking a fixed element g & G 
and considering the space of admissible sequences, which consist of compos- 
able elements gi, ■ ■ ■ ,gN G G such that gi ■ ■ ■ gN = 5- The discrete action 
for an admissible sequence is then taken to be 

N 

Sh{9i,---,9N) = '^Lh (gn), 

n=l 

and discrete Euler-Lagrange equations are obtained by finding sequences 
which are critical for this action function. In the case of the pair groupoid 
Q X Q, fixing an element of the groupoid corresponds simply to fix- 
ing the endpoints g = {qo,qN), while the set of admissible sequences 
(^0) 9i) J • • • J iqN-i,qN) & Q X Q can be identified with the sequence of con- 
figuration points qo, . . . , qj\f £ Q. 

Marrero et al. (2006) showed that the discrete Euler-Lagrange equations 
can be expessed in terms of left- and right-invariant vector fields on G, 
each of which can be identified with sections of the Lie algebroid AG Q 
associated to G. To describe this, we must first introduce some notation. 
Let a,/3: G ^ Q denote the source and target maps, e: Q ^ G denote the 
identity section, and i: G ^ G denote the inversion map on G. For any 
g & G, define the left- and right-translation maps, respectively, by 

Eg : a-' {Pig)) ^ a'^ {a{g)) r, : P'^ {a{g)) ^ P'^ {P{g)) 

g' ^ 99' 9' ^ g'g- 

Then, given a section X eV (AG) , the left-invariant vector field G X{G) 
is defined by ^{g) = T^(^p(^g^^lg{X {P{g))). Similarly, the right-invariant 
vector field ^ G %{G) is given by t{g) = - {T,^a{g)) (.Tg o i)) iX{a{g))). 
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This implies the foUowing relationship between the brackets -J on T (AG) 
and [•, •] on X{G): 

With these definitions, the discrete Euler-Lagrange equations are 

^ [Lh] (gn) = ^ [Lh] (ffn+i) , n = 1, . . . , iV - 1, 
for all sections X G F (AG) . In the pair groupoid case, one sees that 

^ [Lh] {qn-i, qn) = diLh {qn-1, qn) " X (q^) 
^ [Lh] {qn, qn+i) = -doLh {qn, ^n+l) " X (g„) , 

for any X G X{Q), so this formulation agrees with the earlier expression of 
the discrete Euler-Lagrange equations on Q x Q. 

For the symplectic approach, we begin by defining the cotangent groupoid 
T*G ^ A*G, where the base A*G ^ Q is the dual of the Lie algebroid 
AG — >■ Q. This is a symplectic groupoid, with a canonical symplectic struc- 
ture on T*G and Poisson structure on A*G: moreover, the source and tar- 
get maps are anti- Poisson and Poisson, respectively (Coste et al., 1987; 
Marie, 2005). Explicitly, the source map a: T*G A*G and target map 
/3: T*G A*G may be defined such that, taking e T*G, 

a{fi)-X {a{g)) = ^ • l{g), /3 (m) • X (/3(<7)) = ■ 1({g), 

for every section X gT {AG). (The multiplicative structure of T*G will not 
be necessary here, so we omit discussion of it.) 

Now, the discrete Lagrangian Lh - G ^ M generates a Lagrangian sub- 
manifold of the cotangent groupoid, dLh{G) C T*G. If Lh is suitably non- 
degenerate, it follows that this submanifold determines a Poisson automor- 
phism A*G A*G (Weinstein, 1996; Marrero et al, 2006). Specifically, 
define the discrete Legendre transforms ¥^Lh ■ G A*G by 

F'Lh = ao dLh, ¥+Lh = ^ o dL^. 

Therefore, {{¥~ Lh{g) Lh{g)) | g € G} is a coisotropic relation on A*G, 
and if ¥~Lh is a diffeomorphism, then this relation is the graph of the 
discrete flow map 

¥+Lh o {¥^Lh)'^ : A*G A*G. 

To see that this is compatible with the discrete Euler-Lagrange equations, 
observe that for any X eT {AG) , 

{a o dLh) {g) ■ X {a{g)) = dLh{g) ■ J^{g) = 1 [Lh] {g) 
{P o dLh) {g) ■ X {P{g)) = dLh{g) ■ 1c{g) = ^ [Lh] {g). 
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Thus, the discrete Euler-Lagrange equations correspond to the condition 
F+L/i (gn) = F~Lft (5n+i)- Equivalently, this means that gi, . . . ,gN is a solu- 
tion of the discrete Eulcr-Lagrange equations when dL^ (gi) , . . . ,dLh (^at) 
is a composable sequence in T*G. 



For the continuous Hamilton-Pontryagin principle, the key idea was to relax 
the condition for curves in TQ to be second-order, but to enforce this condi- 
tion weakly using Lagrange multipliers in T*Q. Analogously, we will develop 
a discrete Hamilton-Pontryagin principle by relaxing the requirement that 
the sequence gi,--- ,gN £ Gis admissible — requiring only that it be admissi- 
ble up to homotopy — and will weakly enforce the condition that this homo- 
topy is, in fact, constant. To allow for the "Lagrange multipliers," these 
homotopies will be given by paths in T*G rather than in G. 

To simplify the exposition, we will first sketch this formulation, in sec- 
tion 3.1, using cotangent bundle coordinates (5,/x) G T*G. The fully intrin- 
sic treatment on T*G will be given subsequently, in section 3.2, along with 
definitions and the proof of the main theorem. 

3.1. Formulation in cotangent bundle coordinates. Suppose that 
{gn, l^n) '■ [0; 1] — >■ T*G, for n = 1, . . . ,N, is a sequence of paths in the cotan- 
gent groupoid. Furthermore, require that the initial point (fl'n(O), /i„(0)) 
of each path lie in the zero section, so that //n(0) = Og^io), that 
gi{0),. . . , 5Ar(0) be an admissible sequence for a fixed element g E G. 

Given such a sequence of paths, take the discrete Hamilton-Pontryagin 
action to be 



3. The Discrete Hamilton-Pontryagin Principle 




Variations of this action are 



dSh (5i,Mi> • • ■,gN,IJ'N) ■ {6gi,SiJ,i, . . .,SgN,Si2N) 




= / ii^^^h {gn{s)) - /x^s)) • 6gn{s) + (^/x„(s) • g'^{s)] ds 




Now, the integral terms vanish when {gml^n) is a solution of 
IJ.'^{s) = dLh ignis)), g'^is) = 0. 
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The second equation states that gn{s) = Qn is constant, so /in(s) naoves 
vertically along the fiber T*^G, beginning at /in(0) = Og„- We use this fact 
to solve the remaining equation, 

IJ'n{s) = ^n(O) + S dLh (gn) = S dLh {Qn) , 

SO in particular, = dLh (gn)- Additionally, gi, . . . , g^^ is an admissible 

sequence, since gn = gn{^)-, which was assumed to be admissible. Finally, 
restricting to {gmfJ-n) where these equations are satisfied, the variation of 
the action is given by the remaining boundary terms 

N N 

n=l n=l 

However, the restricted variations Sgn arc no longer completely arbitrary, 
since they must be tangent to the space of admissible sequences gi, . . . ,gN. 
Therefore, we are back in the case considered by Marrero et al. (2006), so 
these terms vanish precisely when gn is a solution of the discrete Euler— 
Lagrange equations — or equivalently, when . . . , /Xjv(l) € T*G is a 

composable sequence. 

In summary, dSh = implies that gi{l), . . . , g^i^) is an admissible 
sequence in G, that /Un(l) = dL/i(5„(l)), and that . . . , /iiv(l) is a 

composable sequence in r*G. 

3.2. Intrinsic formulation. We will now make precise the approach 

sketched out in the previous section, doing this intrinsically on T*G rather 
than in cotangent bundle coordinates. Let ttq- T*G G he the cotangent 
bundle projection, 6 be the canonical 1-form on T*G, and Co = —dO be the 
canonical symplectic 2-form. 

Given a Lie groupoid G ^ Q, define the fibration E ^ G consisting of 
paths in T*G beginning at the zero section, 

^ = {7: [0, 1] ^ r*G I 7(0) G Og} , 
where the projection is given by 7 1— )• (ttg o 7) (1). 

Definition 3.1. Given a fixed g & G, & sequence of paths 71, . . . ,77V G E 
is said to be admissible if the initial elements {ttq o 7i) (0), . . . , (ttg o 7iv) (0) 
form an admissible sequence in G. 

This implies that the projected sequence (ttg o 71) (1), . . . , (ttg o 7jv) (1) G 
G is only admissible up to homotopy, which is weaker than the usual assump- 
tion that gi, - ■ ■ ,gN niust actually be an admissible sequence. 
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Definition 3.2. Let L/j : G ^ M be a discrete Lagrangian, and 71, ... , € 
E be an admissible sequence of paths. Then the discrete Hamilton- 
Pontryagin action of this sequence is 



The sequence is said to satisfy the discrete Hamilton-Pontragin principle if 



Theorem 3.3. Let ji, . . . ,ji\f E E be an admissible sequence of paths, and 
denote fin = 7n(l) G T*G and gn = ttg (fJ-n) G G for n = 1, . . . , N . If the 
sequence satisfies the Hamilton-Pontryagin principle, then 

(i) gi, . . . ,g]y is an admissible sequence in G, 
(a) jdn = dLh {gn) forn = l,...,N, 
(Hi) /xi, . . . ,/xjv is a composable sequence in T*G. 

These properties may be thought of as, respectively, (i) the discrete 
second-order curve condition, (ii) the discrete Legendre transform, and (iii) 
the discrete Euler-Lagrange equations. 

Proof. Take a variation of the Hamilton-Pontryagin action to obtain 



Applying Cartan's "magic formula," the Lie derivative terms become 




dS'/i(7i,...,7iv) = 0. 



dS'^(7i,...,7iv) ■ ((57i,...,(57;v) 





Therefore, the variation of the action is 



dS'/i(7i,...,7jv) • ((571, . . . ,(57jv) 




For variations taken along the fibers of E, this is stationary when 
h'n(s)^ = -'^hdLh, n=l,...,N, 
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i.e., the paths 7„ are solutions to Hamilton's equations for the singular 
Hamiltonian H = —L^ o -kq: T*G R. However, since H is constant on 
fibers of T*G (which form a Lagrangian foliation), the Hamiltonian flow is 
fiber-preserving, and hence the projection (ttg ° 7n) (*) = Qn is constant. 
This proves (i), since Qn = [t^g ° In) i*^), which was assumed to be an 
admissible sequence in G. Furthermore, the Hamiltonian vector field is 
constant on fibers of T*G, and it follows that 7n(s) = sdLh {gnj- Therefore, 
/^n = 7n(l) = <^Lh {gn), which proves (ii). 

Finally, restricting the action to these solutions, the integral terms vanish, 
and the remaining boundary terms are 

N N N 

n=l n=l n=l 

This is precisely the variation of the usual discrete action, so these terms 
vanish when the discrete Euler- Lagrange equations are satisfied. Therefore, 
(3 (fin) = a (/Un+i) for n = 1, . . . , — 1, SO ^ui, . . . , jUjv is a composable 
sequence in T*G, which completes the proof of (iii). □ 

Remark 1. For a single time step A'^ = 1, the action 

Sh-.E^R, -f^ f {Lh0 7rG)h{s))ds+ fo 
Jo 

can be understood as a Morse family, which generates the usual Lagrangian 
submanifold dLh{G) C T*G. This is closely related to the work of Milinkovic 
(1999), who studied a similar action principle on T*Q, in connection with 
the Morse homology of generating functions for Lagrangian submanifolds. 
In Milinkovic's formulation, the Lagrangian submanifold is determined by 
the time-1 Hamiltonian isotopy of which in this case, takes the zero 
section to dLh{G). 

4. Example: The Pair Groupoid 

Let G = Q X Q he the pair groupoid, and L/j: Q x Q — >■ M be a discrete 
Lagrangian. Fixing the endpoints g = (goj In)-, suppose that 71, ... , 77V G E 
is an admissible sequence of paths, where 7n(s) = (9n('5))Pn('S)) 'J'n('^))Pn('5)) • 
Since the sequence is admissible, this implies in particular that = 9i(0), 
5n(0) — 1n+ii^) for n = 1,...,A^ — 1, and q\j{^) = Qn- The discrete 
Hamilton-Pontryagin action of this sequence is then 



Sh{ni,---,iN) 

n=l-^0 L 



ds. 
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Taking variations of this action gives 



dSh (71, . . . , 7Ar) ■ ((571, . . . , Sjn) 



n=l 







-pys).6q'J{s)+pi{s).6ql,'{s) 

-6p'iis)-q'J{s) + Spiis)-q},'is)\ds, 
and integrating by parts, this simphfies to 
dS'/j(7i,...,7Ar) • {dji,...,SjN) 



N 



= E 

n=l 



{doL, {q',is),qUs))+p'n'is))-dq',is) 

+ {diLh ((?°(s),gi(s))-pi'(5))-<5gi(6 
- Sp'iis) ■ g°'(s) + Spl,{s) ■ ql'{s) 

p'i{l)-6q^M)+Pi(.l)-Sqi{^) 



ds 



Pnis) 

pi'is) 



For variations taken along the fibers of E, this is stationary when the integral 
terms vanish, so 

-doLh{q^nis),qUs)), Qn'(s) = 0, 
diLh{q%s),qi{s)), ql'{s) = 0, 

for n = I, . . . , N. Therefore, q^is) = q^ and (/^(s) = q^ are constant, so 
PniX) = -dGLh{q^,ql) , pii'f) = diLh {q^^,ql) , n=l,...,N. 

Furthermore, the admissibility assumption implies that qo = qi, qn = 9n+i 
for n = 1, . . . , A'^ — 1, and q]^ = q^- 

Finally, restricting to these solutions, the restricted variations must then 
satisfy Sq^ = 0, Sq}^ = 6q^^^^ = for n = 1, . . . ,N — 1, and dqj^ = 0. Therefore, 
the remaining terms of the action are 

AT AT-l 

5] (-^,0(1) ■ Sq'^+pl,{l) . Sqi) = ■ Sql 

n=l n=l 

which vanish when p\{^) = for n = 1, . . . , A'" — 1. 

In summary, if 71 , ... , 77V satisfies the discrete Hamilton-Pontryagin prin- 
ciple on Q X Q, then the following is true: 

(i) qo = q^, q\ = for n = 1, . . . , AT - 1, and q]^ = qN, 

(ii) (1) = -doLh {ql qi) and (1) = diL^ {q^ q^) forn = 1, . . . , AT, 
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(iii) p^(l)=pO+,(l)forn = l,...,iV-l. 

These are, respectively, the discrete second-order curve condition, the dis- 
crete Legendre transform, and the discrete Euler-Lagrange equations for 
systems on Q x Q. 

Remark 2. The cotangent paths p^, p\ can be seen as Lagrange multiphers, 
which serve to enforce the coniposabihty conditions q\ = q^+i^ 'we\[ as 
the endpoint constraints qi = qo and g]y = q^- Leok and Ohsawa (2008) 
showed that, when Q is a vector space (or given a particular local coordinate 
neighborhood), this can be done even more directly, by simply subtracting 
these points rather than taking a path between them. This leads to an 
alternative choice of the action sum, 

N N-1 
n=l n=l 

whose variations are computed to be 

N 

X [{doLh {qn, qi) + Vn-i) ■ Sq^ + [diLh (g°, qi) - Vn) ■ ^qi] 
n=l 

N-1 

+ Spo ■ {qi - ^o) + X ^P'^ ■ ~ 'iri) + ^PN ■ {qh - qN) ■ 

n=l 

Therefore, this vanishes when 

Pn-i = -doLh {qn, qi) , Pn = diLh {qn, qi) , n=l,...,N 
and when 

qi = qo, qh = qN, qi = qn+i, n=l,...,N -l, 

which are consistent with the results obtained in this section. 

Effectively, the Leok-Ohsawa approach can be interpreted as using a 
smaller fibration E' = Q x Q x T*Q x T*Q, consisting only of the endpoints 
of the paths in to define a Morse family over QxQ. However, when Q is 
not a vector space, there is no intrinsic, global meaning to subtracting two 
points in Q. This action may still be defined locally, in a neighborhood of 
the diagonal of Q x Q, given a choice of local coordinates or a retraction. In 
general, though, one requires the larger space of paths E in order to define 
the discrete Hamilton-Pontryagin principle globally. 
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